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Abstract. This paper considers the problem of modeling the light polarization that 
emerges from an astrophysical plasma composed of atoms whose excitation state is 
significantly influenced by the anisotropy of the incident radiation field. In particular, it 
highlights how radiative transfer simulations in three-dimensional models of the "quiet" 
solar atmosphere may help us to probe its thermal and magnetic structure, from the near 
equilibrium photosphere to the highly non-equilibrium upper chromosphere. The paper 
finishes with predictions concerning the amplitudes and magnetic sensitivities of the 
linear polarization signals produced by scattering processes in two transition region 
lines, which should encourage us to develop UV polarimeters for sounding rockets and 
space telescopes with the aim of opening up a new diagnostic window in astrophysics. 



1. Introduction 

Consider a line of sight (LOS) that makes an angle 6 = 45° with respect to the solar ra- 
dius through the observed point (i.e., \i - cos 6 - 0.707) and let us make the following 
virtual measurements of the emergent Stokes profiles in a "non-LTE line". Assume that 
the observed plasma structure is strongly magnetized (with magnetic strength B - 1200 
G), and consider two different magnetic field orientations: (1) vertical and (2) horizon- 
tal (see Fig. 1). The dotted lines of Fig. 2 show the emergent Stokes profiles produced 
by the Zeeman effect, which for this LOS turn out to be identical for both field orien- 
tations because of the well-known azimuth ambiguity of the Zeeman effect. However, 
the Stokes profiles produced by the Sun are very different. Actually, the Stokes profiles 
corresponding to the vertical field case are given by the solid lines of Fig. 2, while those 
corresponding to the horizontal field case are given by the dashed lines. 




Figure 1. The LOS and the magnetic field orientations considered in Figs. 2 and 3. 
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Figure 2. Emergent Stokes profiles of the He i 10830 A triplet for B = 1200 
G, calculated with the spectral sy nthesis option of the code H AZEL (from HAnle 
and ZEeman Light) developed bv lAsensio Ramos et alJ (l2008h . The "ATOMPOL" 
profiles have been obtained taking into account the joint action of atomic level polar- 
ization and the Hanle and Zeeman effects, while the "NO ATOMPOL" profiles are 
produced by the Zeeman effect alone. The reference direction for positive Stokes 
Q is the perpendicular to the plane formed by the solar radius vector through the 
observed point and the LOS. 



We may thus note that 



• The polarization of the Zeeman effect only depends on the orientation of the 
magnetic field with respect to the LOS. 

• The polarization produced by the Sun depends on both the orientation of the 
magnetic field with respect to the LOS and on its inclination with respect to the 
local vertical. 



Why does the Sun behave this way? Because we are dealing here with an astro- 
physical plasma composed of atoms whose excitation state is significantly influenced 
by the anisotropy of the pumping radiation field. The ensuing radiative transitions pro- 
duce population imbalances and quantum coherences between the magnetic sublevels 
pertaining to the lower and/or upper levels of the spectral line under consideration (i.e., 
atomic level polarization), in such a way that the populations of magnetic substates 
with different values of \M\ are different. This is termed atomic level alignment. The 
radiatively induced polarization of the atomic levels produces selective emission and/or 
selective absorption of polarization components, which give rise to linear polarization 
in the emergent spectral line radiation even in the absence of magnetic fields (e.g., 
Trujillo Bueno 2001, for a detailed review). Via the Hanle effect the atomic level po- 
larization is sensitive to the presence of a magnetic field inclined with respect to the 
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Figure 3. As in Fig. 2, but for B = 100 G. 



symmetry axis of the incident radiation field, so that it is natural to find in Fig. 2 very 
different results for the vertical and horizontal magnetic field cases. 

What kind of Stokes profiles would we observe if the magnetic strength were much 
lower (say, B = 100 G) ? As seen in Fig. 3, if B = 100 G the circular polarization of the 
He i 10830 A triplet is still dominated by the longitudinal Zeeman effect, but the linear 
polarization is fully controlled by atomic level polarization and the Hanle effect with 
no significant influence of the transverse Zeeman effect. Interestingly, this is precisely 
the situation we find in quiet regions of the solar atmosphere for many spectral lines, 
not only for the He i 10830 A triplet considered in Figs. 2 and 3. An observational 
example is shown in Fig. 4, which corresponds to high- sensitivity spectropolarimetric 
measurements of the Ca n IR triplet obtained in collaboration with R. Ramelli (IRSOL) 
using the Zurich Imaging Polarimeter (ZIMPOL) at the THEMIS telescope. This is 
because for weak magnetic fields (with 5^100 G) the Zeeman splitting of the <x and n 
components of many spectral lines is negligible compared with the spectral line width. 

This paper highlights how radiative transfer modeling of the scattering polariza- 
tion observed in the "quiet" solar atmosphere allows us to explore its thermal and mag- 
netic structure, but it does not provide information on numerical methods. Such in- 
formation can be found in several recent reviews. For the case of partial frequency 
redistribution ( PRD), but assuming a two-leve l model atom with no lower-level po- 
larization, see iNagendra & Sampoornal (H()09). For the case of realistic multilevel 
atoms with the possibility of atomic level polariz ation in all the l evels, but assum- 
ing complete frequency redistribution (CRD), see iTrujillo Buenol (l2003h . Unfortu- 
nately, a general PRD theory for mo deling scatterin g line polarization with multi- 
level atoms is not yet available (e. g.. iBelluzzil 1201 1 | ). Nevertheless, the C RD the- 
ory described in the monograph by lLandi Degl'Innocenti & Landolfil (120041 ) is suit- 
able for modeling the polarization observed in the Sun's continuous radiation (see 
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Figure 4. An example of our recent spectropolarimetric observations of the Ca n 
IR triplet in a quiet region close to the solar limb, using ZIMPOL at the telescope 
THEMIS of the Observatorio del Teide. Note that the sign of the U /I signals fluc- 
tuates along the spatial direction of the spectrograph's slit, which suggests spatial 
variations in the azimuth of the magnetic field in the quiet solar chromosphere. The 
reference direction for Stokes Q is as in Fig. 2. From Trujillo Bueno et al. (201 1). 



Section 2.1) and in many diagno stically important spectral lines, suc h as the He i 



10830 A and D 3 multiplets (e.g.. ICasini & Landi Pegl'Innocenti! 120071). the Sr i line 



at 46 07 A (see Sections 2.2 and 2.3), the lines of ionized cerium (|Manso Sainz et al 



2006).Jie_CN lines (e.g.. lAsensio Ramos & Truiillo Bu eno 2003), the Li i 6708 A line 



(jBelluzzi et alJl2009h. the Mg i Min e s dTruiillo B ue noll200ll), the IR triplet of Ca n 
dManso Sainz & Trujillo Buenoll2003l . l2010h and H„ (IStepan & Trujillo Buenoll2010h . 
Moreover, it can also be used for estimating the line-center scattering polarization am- 
plitudes expected in strong resonance lines for which frequency correlations between 
the incoming and outgoing photons are indeed significant (e.g., the Ly^ line investiga- 
tion of Section 3.1). However, PRD effects have to be taken into account in order to 
model the shapes of the linear po larization profiles produced by scattering in strong res- 
onance lines (e.g jAnusha et aljfeoill ). Fortunately, the lower level of several resonance 
line transitions cannot carry atomic alignment because its total angular momentum is 
// = or // = 1 /2, so that PRD modeling efforts based on the standard two-level atom 
approach are of interest (e.g., the Mg n k line investigation considered in Section 3.2). 



2. 3D modeling of scattering polarization in the quiet photosphere 

This section considers some radiative transfer investigations of scattering polarization 
in the following three-dimensional (3D) models of the quiet solar photosphere: (1) a 
sn apshot model taken fro m the hydrodynamical simulations of solar surface convection 
by lAsplund et al.1 (120001 ) (hereafter, the HD model) and (2) a snapshot model taken 
from the magnetoconvection simulations with surface dynamo action carried out by 
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Vogler & SchussleJ (120071) (hereafter, the MHD model). The main aim here is to review 



how the modeling of observations of the intensity and linear polarization of the Sun's 
continuous radiation and of the Sr 1 4607 A line allows us to probe the thermal structure 
and the magnetization of the quiet solar photosphere. 

In both radiative transfer problems "zero-field" dichroism (see Ma nso Sainz & Trujillo Buenol 



2003) can be neglected, so that the 4x4 propagation matrix of the Stokes- vector trans- 



fer equation is diagonafl. Therefore, the relevant equation for the Stokes parameter X 
(with X = I,Q, U) at a given frequency v and direction of propagation Q is given by 

^-X = X - S x , (1) 
or 

where dr = -r/jds defines the monochromatic optical distance (r) along the ray, with s 
being the corresponding geometrical distance and r\i the absorption coefficient. In Eq. 
(1) the expressions of the source-function components are 

S 7 = rS 7 * + (l-r)£ v , (2) 
S Q = rS* Q , (3) 

Su = rS* U9 (4) 

where r is a problem-dependent opacity ratio to be given below and B y the Planck 
function. Moreover, 

S 7 * = S° Q + wfj^Of - l)Sl - V3/i VW?(cos^Re[5?] - sin^Im^]) (5 ) 
+^(1 -fi 2 )(cos2 X Ro[S 2 2 ] - sin2xIm[S^])}, 

S Q=- w Zi{^^ 2 ~ V s o ~ ^ VW?(cos^Re[5?] - sin^Im^]) (6) 
-%\ + y u 2 )(cos2*Re[^] _ sin2^Im[52])), 

and 



S \j = - wf uJi V3( VW?(sin*Re[S ?] + cos*Im[S \]) (7) 
+Msin2^Re[52] + cos2^Im[5 2 ])), 

where 6 = arccos (p) and x are the inclination with respect to the solar local vertical 
and azimuth of the ray, respectively, wf uJ( is a coefficient that is equal to unity for 



! For the spectral line case the propagation matrix is diagonal only if, in addition, the Zeeman splitting of 
cr and 7i components is negligible compared with the spectral line width. This is a good approximation for 
modeling the linear polarization amplitudes of the Sr 1 4607 A line observed in the quiet Sun. 
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the problems considered in the following subsections, and the equations that govern 
the quantities S 2 , Re[S 2 ], Im[S 2 ], Re[S 2 ] and Im[S 2 ] will be given below for 
each particular problem considered. We point out that in these equations the reference 
direction for Stokes Q is the perpendicular to the plane formed by the ray's propagation 
direction and the vertical Z-axis. 



2.1. Polarization of the Sun's continuous spectrum 

The dominant contribution to the linear polarization of the visible continuous spectrum 
comes from Rayleigh scattering by neutral hydrog en in its ground state (Lyman scat- 
tering) and Thomson scattering by free electrons (IChandrasekhad 1 1 960h Fl. The con- 
tribution of these processes to the total absorption coefficient (77/) is quantified by 
cr c = crjNe + o-Rni(H), where crj - 6.653xl0" 25 cm 2 is the Thomson scattering 
cross section, N e the electron number density, n\{H) the population of the ground level 
of hydrogen, and ctr the wavelength-dependent Rayleigh cross section. The contin- 
uum absorption coefficient is given by r\i - k c + <x c , where k c contains all the relevant 
non- scattering contributions (e.g., the bound-free transitions in the H" ion). The expres- 
sions of the Sx source function components can be found easily by taking r = 

and S° = S 2 = / 2 , Re[S 2 ] = Re[ / 2 ], Im[S 2 ] = -Im[/ 2 ], Re[S 2 ] = Re[/ 2 ] and 
Im[S 2 ] = -lm[jj] inEqs. (5)-(7) (e.g. JTrujillo Bueno & Shchukina 2009). 

Note that the Jq quantities (with K = 0, 2 and Q = 0, 1,2) are the spherical 
comp onents of the radiation field tensor (see §5.11 in lLandi Pegl'Innocenti & Landolfil 



2004), which quantify the symmetry properties of the radiation field at the spatial point 

under consideration. Thus, J® is the familiar mean intensity, j> ^^=(3// 2 - l)/ y ^ 

quantifies its anisotropy, while the real and imaginary parts of J 2 Q (with Q = 1,2) 

(i.e., Re[/g] and Im[/^], respectively) measure the breaking of the axial symmetry. 

Obviously, the real and imaginary parts of j\ and / 2 are zero in a plane-parallel or 
spherically symmetric model atmosphere, but they can have significant positive and 
negative values in a 3D stellar atmospheric model (see figure 3 of Trujillo Bueno & 
Shchukina 2009). Therefore, Eqs. (6) and (7) tell us that U - in ID models, and 
that the key observational signatures of the symmetry breaking effects in a 3D model 
(caused by its horizontal atmospheric inhomogeneities) are non-zero U signals at any 
on-disk position and non-zero Stokes Q and U signals for the LOS with fi = 1, which 
corresponds to forward- scattering geometry. 

Figure 5 shows the Q/I and U/I images that we would see at 4600 A if we could 
observe with high polarimetric precision the solar continuum polarization at very high 
spatial and temporal resolution. Obviously, without spatial and/or temporal resolution 
and the only observable quantity would be Q/I, whose wavelength variat ion at 
a solar disk position close to the limb has been determined semi-empirically by lStenflol 
(l2005h . Note that the spatially-averaged Q/I of the solar continuum radiation at each 
\i position is largely determined by the effective polarizability, r - <t c /(k c + <x c ), and 
by the fractional radiation field's anisotropy, J^/Jq, which in turn depend on the den- 
sity and thermal structure of the stellar atmosphere under consideration. Therefore, the 



2 The H~ ion cannot produce polarization through Rayleigh scattering because it has a single bound state, 
Is 2 ^o. Moreover, there is no contribution from "zero-field" dichroism because the H~ level cannot be 
polarized. 
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Figure 5. The continuum scattering polarization signals at 4600 A calculated for 
three line of sights in Asplund et al.'s (2000) hydrodynamical model of the solar 
photosphere (i.e., the HD model), accounting for the diffraction limit effect of a 1-m 
telescope. Although the polarization amplitudes are very low at 4600 A, they are 
one order of magnitude larger at 3000 A. The solid line contours in the right panel s 
delineate the upflowing granular regions. From lTrujillo Bueno & Shchukind (12009). 



more realistic the thermal and density structure of a given solar atmospheric model is, 
the closer to the empirical data will be the calculated linear polarization of the s olar con- 
tinuum radiation. Interestingly, Fig. 12 of iTrujillo Bueno & Shchukinal (l2009h demon- 
strates that 3D radiative transfer modeling of the polarization of the Sun's continuous 
spectrum in the above-mentioned 3D hydrodynamical model of the solar photosphere 
shows a notable agreement with Stenflo's (2005) semi-empirical data, significantly bet- 
ter than that obtained via the use of ID semi-empirical models. 

2.2. The polarization of the Sr i 4607 A line and the Hanle effect of a microtur- 
bulent field in a 3D hydrodynamical model of the quiet photosphere 

In this and in the following subsection we focus on the diag nostic problem of the mag- 
netism of the quiet solar photosphere (e.g., the review by ISanchez Almeida! 120041 ). To 
this end, we need to measure and interpret polarization in spectral lines. It is important 
to note that determining the mean magnetization of the quiet Sun requires finding how 
much flux resides at small scales and that, to this end, it is crucial to measure and inter- 
pret the linear polarization produced by atomic level polarization and its modification 
by the Hanle effect. 

If at the spatial grid point under consideration the magnetic field of strength B is 
assumed to have an isotropic distribution of orientations within a volume X 3 , with £ 
the mean-free-p ath of the line photons, the equations that govern the Sq quantities of 
Eqs. (5)-(7) are dTmjillo Bueno & Manso Sainzlll999h : 
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and 



S° = (1-6)7° + eB v , 



(8) 



/ s 2 



Re[S 2 ] 
Im[Sf] 



ImtS^] 



W 



(1-6) 



JuJl I + 5 (2) (1 _ 6 ) 



J 2 ^ 



Re[J?] 
-Im[^] 



Re[J5] 



-Im[^] 



(9) 



where 6 can be interpreted as the probability that an atomic excitation event is due to 
inelastic collisions with electrons, 6^ is the upper-level depolarizing rate (in units of the 
Einstein A u \ coefficient) due to elastic collisions with neutral hydrogen atoms, and < H (2) 
is the Hanle depolarization factor whose value is = 1 for B = G and Hrft^ = 1/5 
for B > 2?satur (where B satm is the Hanle saturation field, which is at least 200 G for 
the Sr i 4607 A line). It is important to note that for the spectral line cases considered 
here and in the following sections the radiation field tensors are given by angular and 
frequency-weighted averages of the Stokes parameters, with the Voigt profile y as the 
frequency-dependent weighting function. For this reason, in Eqs. (8) and (9) they are 
instead denoted with the symbol Jq. It is also worthwhile pointing out that for the 

spectral line case S q = -^==p^, with being the multipolar components of the 

upper-level density matrix normalized to the overall population of the transition's lower 
level. Therefore, Eq. (9) describes the transfer to the atomic system of the symmetry 
properties of the radiation field, and note that it is the largest for line transitions with 
wf > Ji - 1 and when in the absence of magnetic and collisional depolarization. For the 

Sr i 4607 A line // = 0, J u = 1 and = 1. Moreover, in Eqs. (2)-(4) we have now 
r = Ki(f) y /(Ki(f)y + k c ), with ki the line-integrated opacity. 

Figure 6 shows the center-to-limb variation of the Q/I and U/I line-center signals 
of the Sr i 4607 A line, which we have obtained by solving the scattering line polar- 
ization problem in the above-mentioned HD model, without including any magnetic 
field, but accounting for the diffraction limit effect of a lm telescope. Note that the 
calculated linear polarization amplitudes have sizable values and fluctuations, even at 
the very center of the solar disk, where we observe the forward scattering case (e.g., the 
standard deviations of the Q/I and U /I fluctuations are approximately 0.3% atyu = 0.5). 
Trujillo Bueno & Shchukinal (l2007h pointed out that the predicted small-scale patterns 
in Q/I and U/I are of great diagnostic value, because they are sensitive to the thermal, 
dynamic and magnetic structure of the quiet solar atmosphere. Therefore, it would be 
important to try to observe them with an imaging polarimeter attached to a suitable 
telescope, such as that of the SUNRISE balloon-borne telescope. 

For the moment, the scattering polarization observations of the Sr i 4607 A line 
that have been published lack spatial resolution (e.g., see the data points of Fig. 7). 
Interestingly, when the calculated /, Q and U profiles at each jjl solar-disk position are 
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Figure 6. The fractional linear polarization amplitudes of the Sr i 4607 A line 
calculated in the HD model of the solar photosphere, neglecting any interaction with 
the continuum polarization because at /14607 it is almost two orders of magnitude 
smaller, as shown in Fig. 5. From Trujillo Bueno & Shchukina (2007). 



spatially averaged over scales larger than those of the solar granulation pattern we ob- 
tain U/I^O and the Q/I vs. \i line-center amplitudes given by the green dotted line of 
Fig. 7. The very significant di screpancy between the obse rved and the calculated scat- 
tering polarization signals led lTrujillo Bueno et al.l (l2004h to the following conclusion: 
the internetwork regions of the quiet solar photosphere are permeated by a spatially un- 
resolved tangled magnetic field, whose mean field strength at heights /?~3()0 km above 
the visible solar surface is (2?)^ 130 G when no distinction is made between granular 
and intergranular region^. This first conclusion was reached by modeling the ob- 
served scattering line polarization through the self-consistent solution of Eqs. (8), (9) 
and (1) in the above-mentioned hydrodynamical model. Note that such a 3D photo- 
spheric model is unmagnetized and that, therefore, our (5)^130 G conclusion is based 
on the following approximations for the description of the quiet Sun magnetic field that 
produces Hanle depolarization: (1) that the magnetic field is tangled at subresolution 
scales, with an isotropic distribution of directions and (2) that the shape of the proba- 
bility distribution function, describing the fraction of quiet Sun occupied by magnetic 
fields of strength B, is exponential (i.e., PDF(£) = ^~ BI{B) ). 



3 Note that (5)^130 G is an order of magnitude larger than the ol d Hanle-effect estimations b y Stenflo 
(1982) and Faurobert-Scholl et al. (1995). For a detailed review see lTruiillo Bueno et alJ (l2006h . 

4 It is of interest to men tion that in Cattaneo's (1999) numerical experiments of surface dynamos 
Thelen & Cattaneo (2000) noted how cleanly the PDF track an exponential behavior once the boundary 
effects are removed. 
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2.3. The polarization of the Sr 1 4607 A line and the Hanle effect in a 3D magneto- 
convection model of the quiet photosphere with surface dynamo action 

If at the spatial grid point under consideration the magnetic field has a given strength 
B, inclination 6b and azimuth xb (e.g., as is the case with the magnetic field of magne- 
toconvection simulations) the equation s that govern the S q quantities of Eq. (5)-(7) are 
dManso Sainz & Truiillo Buendll999h : 



+ eB v 



(10) 
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where T u = 8.79 x 10 6 g u B/A u i (with g M as the upper-level Lande factor, B as the mag- 
netic strength in gauss and A u \ as the Einstein coefficient for spontaneous emission in 
s" 1 ). The -coefficients of the magnetic kernel M depend on the inclination (6 B ) 
of the local magnetic field vector with respect to the vertical Z-axis and on its azimuth 
(Xb)> These equations for the S q multipolar components have a clear physical meaning. 
In particular, note that the magnetic operator M couples locally the K = 2 components 
among them (the Hanle effect). 

The black dotted line of Fig. 7 shows the s patially averaged scattering po l arizat ion 
amplitudes of the Sr 1 4607 A line calcul a ted bylShchukina & Trujillo Buenol (l201lh in 
the MHD model of IVogler & Schussled (12007b . but after imposing B = G at each 
spatial grid point. Interestingly, the resulting scattering polarization amplitudes are 
similar (but not identical) to those computed in the HD model (see the green dotted 
line). The differences between the two curves are due to the fact that the thermal and 
density structure of the two snapsho t models are not fully identical. The black dashed 
line shows the Q/I amplitudes that IShchukina & Trujillo Buenol (|201lb obtain when 
taking into account the Hanle depolarization produced by the actual magnetic field of 
the MHD surface dynamo model, whose magnetic Reynolds number is /? m ^2600. The 
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Figure 7. Center-to-limb variation of the (spatially-averaged) Q/I scattering am- 
plitudes of the photospheric line of Sr i at 4607 A. The data correspond to various 
observations taken by several authors during a minimum and a maximum of the 
solar activity cycle. The two green lines sho w scattering polarizat ion amplitudes cal- 
culated in the 3D hydrodynamical model of Asplu ncTet al.l (|2000), without including 
any magnetic field (green dotted line) and including the Hanle depolarization of a 
microturbulent field with an exponential PDF characterized by a mean field strength 
(#)=130G (green dashed line). The two b lack lines show scat tering polarization 
amplitudes calculated in the MHD model of Vogler & Schussler (|2007|) . neglecting 
its magnetic field (black dotted line) and taking into account the Hanle depolariza- 
tion of its magnetic field (black dashed line). We point out that the observations 
can be approximately fitted by multiplying each grid-point magnetic strength by a 
scaling factor, such that at heights /z^300 km the resulting (B) is approximately 130 
G. The red solid line shows the calculated scattering polarization amplitudes when 
imposing B > ff S atur at each g rid point in the MHD model. For more information see 
IShchukina & Trujillo Buenol (1201 ll) . 



fact that the resulting Hanle depolarization is too small to explain the observations is 
not surprising because the mean field strength at a height of about 300 km in the MHD 
snapshot model is onl}0 (B)&15 G (i.e., about an order of magnitude smaller than 
the (B) ^130 G value inferred by Trujillo BuenoetaL In order to investigate 

whether we can fit the observations of Fig. 7 with a magnetic field topology identical 
to that of the MHD surface dynamo snapshot we have increased only the magnetic 
strength of each spatial grid point by multiplying it by a scaling factor. Remarkably, 
the mean field strengt h that we have to reach this way to fit the observations of Fig. 
7 is about 130 G (see Shchul rina & Trujillo Buenol 1201 lh . Note that this conclusion 



5 Note also that in the idealized local dynamo experiments of Cattaneo ( 1999) the average unsigned field 
at the surface is between 10 and 20 G, if an equipartition strength of 400 G is assumed. 
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Figure 8. Theoretical estimate of the fractional linear polarization signals of the 
Ly a line produced by forward scattering processes in the solar transition region, 
taking into account the Hanle effect caused by the presence of a horizontal mag- 
netic field with the indicated strength (in G). The reference direction for positive 
Stokes Q is the parallel to the magnetic field vector. Note how in forward scat- 
tering geometry (i.e., for a LOS with [i = 1) the fractional linear polarization am- 
plitude of Ly a increases as the magnetic strength of the inclined field is increased. 
The semi-empirica l model u sed for the radiative transfe r calculations is FAL-C (see 
iFontenla et al.ll 19931) . From lTrujillo Bueno et al.l (|201 1|) . 



has now been reached without assuming that the shape of the PDF is exponential and 
without using the microturbulent field approximation. Clearly, both approximations 
were suitable to find the mean field strength of the "hidden" field (Truji llo Bueno et al 



2004). In particular, it is interesting to note that the fact that the red line of Fig. 7 
practically coincides with 1/5 of the zero-field linear polarization amplitudes indicates 
that the microturbulent field approximation is indeed a suitable one. 

Finally, it may be of interest to mention the second conclusion of Trujill o Bueno et al. 



( 20041 ). which resulted from their joint analysis of t he Hanle effect in the Sr i 4607 
A line and in C2 lines of the Swan system (see also iTrujillo Bueno et al.1 12006): the 



downward-moving intergranular lane plasma is pervaded by relatively strong tangled 
fields at subresolution scales, with (B) >2? satur (where 2? satur ~ 200 G is a lower limit to 
the Hanle saturation field of the Sr 1 4607 A line). This conclusion, which implies that 
most of the flux and magnetic energy reside on still unresolved sc ales, is also supported 
by the recent investigation of IShchukina & Trujillo Buenol (1201 ll) . 



3. Scattering polarization and the Hanle effect in transition region lines 

A recent review on spectropolarim etric investigations of the magnetization of the quiet 
Sun chromosphere can be found in lTrujillo Buenol (120101) . In the remaining part of this 
section I focus instead on advancing some results of two recent theoretical investiga- 
tions about the linear polarization signals expected for the hydrogen Ly^ line at 1215 A 
and the Mg 11 k-line at 2795 A, produced by scattering processes in the solar transition 
region plasma itself. 
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3.1. Ly a in the solar transition region 

The intensity profile of the hydrogen hy a line has been measured on the solar disk by 
several instruments on board rockets or space-based telescopes (e.g., by the SUMER 
spectrometer on SOHO). These measurements show that Ly a is always in emission at 
all on-disk positions and times, and that the emission originates in the upper chromo- 
sphere and transition region. The center-to-limb variatio n of the observed Ly a intensity 
profiles in the quiet Sun is very small or negligible (e.g. JCurdt et al.ll2008h . This implies 



that the outgoing radiation shows no significant variation with the heliocentric angle 6. 
However, this does not mean that the illumination of the transition region atoms is 
isotropic - which would imply that scattering polarization in Ly a is zero. The incom- 
ing radiation (i.e., the one which illuminates the transition region atoms from "above") 
calculated by J. Stepan (IAC) and I in semi-empirical models of the solar atmosphere 
shows significant limb brightening, so that atoms in the transition region are indeed 
irradiated by an anisotropic radiation field. The ensuing radiation pumping induces 
population imbalances among the sublevels of the upper level 2p 2 P^/2 of the Ly^ line, 
which is the only hydrogen level that makes a significant contribution to its scattering 
polarization. The atomic polarization of the 2p 2 F?,/2 level produces linear polarization 
in Ly a , which is modified by the presence of a magnetic field through the Hanle ef- 
fect. The Q/I amplitudes we have calculated in various models of the solar atmosphere 
vary between 0.5% and 0.05%, approximately, depending on the strength of the mag- 
netic field and the LOS. For example, Fig. 8 shows the calculated linear polarization 
signals caused by the Hanle effect in the forward- scattering geometry of a disk-center 
observation. Note that in order to be able to detect the presence of a horizontal mag- 
netic field of 30 G in the bulk of the solar transition region the polarimetric sensitivity 
of t he measurement should be at least 0.1%. Interestingly, this Ly^ investigation (see 
also lTrujillo Bueno et al.ll2005h has motivated a proposal to NASA (led by Japan, USA 
and Spain) for the development of a Chromospheric Lyma n- Alpha Spectro-Polarimeter 
(CLASP) for a sounding rocket (see llshikawa et al.ll201lh . 



3.2. Mg ii k in the solar transition region 

The CRD theory on which the radiative transfer modeling of the previous sections is 
based treats the scattering line polarization phenomenon as the temporal succession of 
lst-order absorption and re-emission processes, interpreted as statistically independent 
events. As mentioned in the introduction, this theory can also be used for estimating 
the expected scattering polarization amplitudes in strong resonance lines for which fre- 
quency correlations between the incoming and outgoing photons are indeed significant 
(e.g., see Fig. 8 for the case of hy a in forward scattering geometry). However, it is im- 
portant to emphasize that PRD effects have to be taken into account in order to be able 
to provide reliable predictions on the shapes of the linear polarization profiles produced 
by sca ttering in strong resonance lines like Ly^ and Mg n k. Recently. lSampoorna et al.1 
(l2010h have studied in great detail the sensitivity of PRD scattering polarization profiles 
to various atmospheric parameters. To this end, the authors applied the ve ry efficient ra- 



diative transfer method described in lSampoorna & Trujillo Buenol (l2010h . An example 



of the type of investigations that can be carried out with a PRD radiative transfer code 
based on that numerical method is considered in Fig. 9, which shows the complicated 
shape of the Q/I profile of the Mg n k-line caused by scattering processes with PRD 
effects in a semi-empirical model of the solar chromosphere and transition region. 
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Figure 9. The close to the limb (ja = 0. 1) Q/I profile of the Mg n k-line produced 
by line scattering processes in the FAL-C semi-empiri cal model of the solar atmo- 
spher e, taking into account PRD effects as explained in lSampoorna & T rujillo Buenol 
( 2010) and without accounting for the continuum polarization contribution. The right 
panel shows better the line-core spectral region, indicating the line-center Hanle de- 
polarization produced by a random- azimuth horizontal magnetic field with the indi- 
cated strength. The dotted-line shows the only part of the PRD Q/I profile that can 
be estimated via the CRD theory. The reference direction for Stokes Q is as in Fig. 
2. This figure stems from a collaboration between M. Sampoorna (IIA), J. Stepan 
(IAC) and J. Trujillo Bueno (IAC). 



Off-limb observations of spicules in Ly a and Mg n k at various heights above the 
visible solar limb are expected to show higher polarization amplitudes. Both types of 
observations (on-disk and off-limb) would provide precious information on the physi- 
cal conditions of the solar transition region from the chromosphere to the 10 6 K solar 
coronal plasma, even with a spatial resolution of only 1 or 2 arcseconds. 



4. Concluding comment 

Radiative transfer modeling of solar light polarization is the essential link between the- 
ory and observations. In general, realistic modeling requires solving a very complicated 
and still unsolved problem: non-LTE, multilevel, PRD radiative transfer, taking into ac- 
count the joint action of the Hanle and Zeeman effects in 3D atmospheric models of the 
extended solar atmosphere. Fortunately, the Sun is subtle but not malicious: it allows us 
to learn how it works in a stepwise manner, while continuously motivating us to open up 
new research windows in astrophysics through the development of increasingly novel 
theories, instruments, modeling and data analysis techniques. 
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